UNCLASSIFIED 


AD  NUMBER 

ADB004699 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
only;  Test  and  Evaluation;  MAY  1975.  Other 
requests  shall  be  referred  to  U.S.  Army 
Ballistic  Research  Laboratories,  Attn:  AMXBR- 
SS,  Aberdeen  Proving  Ground,  MD  21005. 


AUTHORITY 

USAARDC  ltr,  Apr  1978 


THIS  PAGE  IS  UNCLASSIFIED 


THIS  REPORT  HAS  BEEN  DELIMITED 
AND  CLEARED  FOR  PUBLIC  RELEASE 
UNDER  DOD  DIRECTIVE  5200,20  AND 
NO  RESTRICTIONS  ARE  IMPOSED  UPON 
ITS  USE  AND  DISCLOSURE, 

DISTRIBUTION  STATEMENT  A 

APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED, 


004699 


REPORT  NO.  1781 

PROPAGATION  CHARACTERISTICS  OF 
NARROW  X-RAY  PULSES 


Heber  D.  Jones 
Donald  Eccleshall 
Judith  K.  Temperley 


May  1975 


Distribution  limited  to  US  Government  agencies  only;  Test  and 

^■A’'  /lj  0ther  re<luests  for  this  document  must  be 

ATTN^  AMyRfiD«eC^r^USA  ®a,1!stic  Research  Laboratories, 
Alin.  AnXBR-SS,  Aberdeen  Proving  Ground,  Maryland  21005. 


USA  BALLISTIC  RESEARCH  LABORATORIES 

ABERDEEN  PROVING  GROUND,  MARYLAND 


ran>  Tv'™™} MS  ■ 

s " $ |:S'wa i . 

m «® 

< *■  \\  m w 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


Secondary  distribution  of  this  report  by  originating 
or  sponsoring  activity  is  prohibited. 

Additional  copies  of  this  report  may  be  obtained 
from  the  Defense  Documentation  Center,  Cameron 
Station,  Alexandria,  Virginia  22314. 


The  findings  in  this  report  are  not  to  be  construed  as 
an  official  Department  of  the  Army  position,  unless 
so  designated  by  other  authorized  documents. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whtn  Dmtm  Rntorod) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

1 REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

REPORT  NO.  1781 

3.  RECIPIENT'S  CATALOG  NUMBER 

A.  TITLE  (mnd  Subttllo) 

Propagation  Characteristics  of  Narrow  X-Ray 
Pulses 

S.  TYPE  OF  REPORT  ft  PERIOD  COVERED 

X July  73  - 30  June  74 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7 AUTHORO) 

Ueber  D.  Jones,  Donald  Eccleshall,  and 
Judith  K.  Temper ley 

8.  CONTRACT  OR  GRANT  NUMBERfn) 

9 PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

USA  Ballistics  Research  Laboratories 
Aberdeen  Proving  Ground,  Maryland  21005 

10.  PROGRAM  ELEMENT,  PROJECT.  TASK 
nrA»£A  * WORK  UNIT  NUMBERS 

RDT5I-  Project  No. 

1T161101 A91A 

U CONTROLLING  OFFICE  NAME  AND  ADDRESS 

US  Army  Materiel  Command 
5001  Eisenhower  Avenue 
Alexandria,  Virginia  22304 

12  REPORT  DATE 

MAY  1975 

13.  NUMBER  DF  PAGES 

34 

U MONITORING  AGENCY  NAME  a ADDRESS^//  dllltront  from  Controlling  Ofllc •) 

IS.  SECURITY  CLASS,  (ol  thlo  report) 

IJnclassi  fied 

iSa.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 

16.  distribution  statement  (oi  im.  Report) 

Distribution  limited  to  US  Government  agencies  only;  Test  and  Evaluation; 
May  1975.  Other  requests  for  this  document  must  be  referred  to  Director, 
USA  Ballistic  Research  Laboratories,  ATTN:  AMXBR-SS,  Aberdeen  Proving 

Ground,  Maryland  21005 

17.  DISTRIBUTION  STATEMENT  (or  thm  mbotr* cl  onlorod  In  Block  70,  II  dtltercnt  from  Report) 

18.  SUPPLEMENTARY  NOTES 

19  KEY  WORDS  fConfinue  on  /•(•*«*•«  •)(!•  11  n»c»*»mry  mid  Identify  by  block  numbor) 

X-Ray 

Hydrodynamics 

Laser 

Propagation 

20.  ABSTRACT  (Continue  on  rmvoreo  eld*  II  n*c***mry  mnd  Identity  by  block  numtxrj 

The  disturbance  produced  in  a gaseous  medium  by  the  absorption  of  an 
intense,  narrow  pulse  of  x-rays  is  such  that  the  absorption  of  later  x-rays 
which  follow  the  same  path  is  greatly  reduced.  Various  processes  which 
affect  this  absorption  reduction  have  been  investigated  and  a numerical 
technique  has  been  employed  to  obtain  quantitative  predictions  of  the 

range.  1 

DD,J^“»  1473  eoitiom  of  i nov  is  obsolete  UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWh#n  Dot*  Bnt*r*d) 


TABLE  OF  CONTENTS 


LIST  OF  ILLUSTRATIONS  

I . INTRODUCTION 

II.  ELEMENTARY  CONSIDERATIONS  .... 
HI.  NUMERICAL  TECHNIQUE 

IV.  MATHEMATICAL  FORMULATION  .... 

A.  Photoelectric  Effect  .... 

B.  Energetic  Electron  Energy  Loss 

C.  Compton  Effect  

D.  Mass  Motion  

E.  Other  Processes  ....... 

V.  PRELIMINARY  RESULTS  

VI.  SUMMARY  

REFERENCES  

DISTRIBUTION  LIST  , 


Page 


5 

7 

8 

10 

11 

11 

14 

16 

17 

19 

22 

27 

29 

33 


3 


lg  — 

H PRECEDXN3  PAGE  BLANK-NOT  FILMED 


"■'WjuV-A'-'w 


1 


LIST  OF  ILLUSTRATIONS 


Figure 


Graphic  representation  of  the  problem.  Part  (a)  shows 
a longitudinal  (temporal)  profile  of  an  x-ray  pulse 
and  part  (b)  shows  a radial  profile.  The  pulse  is 
then  a long,  thin  sliver  and  we  consider  its  propa- 
gation successively  through  short  intervals  as  shown 
in  part  (c) 

X-ray  pulse  shape  and  some  resulting  ion  densities 
within  the  beam  as  a function  of  time  for  the  initial 
propagation  interval  


Remaining  pulse  energy  versus  propagation  distance  in 
the  effusion  approximation  


Temperature  versus  time  in  the  effusion  approximation 
within  the  beam  (region  M)  and  surrounding  the  beam 
(region  W) 


5.  Radial  mass  density  (relative  to  STP  air)  distribution 
at  two  times  for  the  laser  parameters  given  in 
Table  I.  Times  t^  and  t0  are  80  ns  and  60  ns,  re- 
spectively, before  the  arrival  of  the  x-ray  pulse 
peak.  The  radial  x-ray  pulse  shape  is  uniform  and 
extends  to  r = 5 y 


6.  Same  type  of  plot  as  Figure  5 except  that  the  radial 
pulse  shape  is  Gaussian  with  a full  width  at  half 
maximum  of  10  p * 
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I.  INTRODUCTION 


Since  its  birth  less  than  twenty  years  ago,  modern  laser  technology 
has  been  following  a natural  evolutionary  path  to  higher  powers  and 
shorter  wavelengths.  Single  laser  pulses  with  energies  of  hundreds  of 
joules  arc  now  available  and  laser  action  has  been  observed  in  the 
vacuum  ultra-violet  region  of  the  electromagnetic  spectrum3.  Efforts 
arc  also  being  made  to  observe  lasing  of  x-raysl+_‘'.  The  characteristics 
of  possible  x-ray  lasers  are  discussed  in  Reference  10  which  concludes 
that  early  devices  will  have  diameters  of  a few  microns,  short  pulse 
lengths  (10  1 to  10"  13  seconds),  and  low  powers.  More  recently,  means 

to  achieve  longer  pulses  have  been  identified11. 


1.  T.  F.  Stratton,  "A  High  Energy,  Short  Pulse  CO2  Laser  System", 

VIII  International  Quantum  Electronics  Conference,  San  Francisco,  CA, 

June,  1974,  Abstract  B.2. 

2.  A.  Crocker,  H.  M.  Lamberton,  E.  W.  Parcell  and  H.  C.  Probyn, 

500  J CO2  Electron  Beam  Laser  Controlled  by  a Glow  Discharge  Electron 

Gun  , VIII  International  Quantum  Electronics  Conference,  San  Francisco,  CA, 

June,  1974,  Abstract  B.4.  ’ 

3.  A.  Wayne  Johnson  and  J.  B.  Gerardo,  "The  Vacuum-Ultra-Violet 
Xenon  Laser",  VIII  International  Quantum  Electronics  Conference, 

San  Francisco,  CA,  June,  1974,  Abstract  M.2. 

4.  M.  A.  Duguay  and  P.  M.  Rentzepis,  "Some  Approaches  to  Vacuum 
UV  and  X-Ray  Lasers",  Appl.  Phys.  Letters,  Vol.  10,  350  (1967). 

5.  R.  C.  Elton,  R.  W.  Waynant,  R.  A.  Andrews,  and  M.  H.  Reilly, 

"X-Ray  and  Vacuum-UV  Lasers",  NRL  Report  7412,  May,  1972. 

6.  P.  J.  Mallozzi,  et  al.,  "X-Ray  Emission  from  Laser-Generated 
Plasmas",  ARPA  Technical  Report  1723,  January,  1972. 

> 

7.  M.  0.  Scully,  W.  II.  Louisell,  and  W.  B.  McKnight,  "A  Soft 

X-ray  Laser  Utilizing  Charge  Exchange",  VIII  International  Quantum 

Electronics  Conference,  San  Francisco,  CA,  June,  1974,  Abstract  M.9. 

8.  R.  A.  McCorkle  and  J.  M.  Loyce,  "Threshold  Conditions  for 

Amplified  Spontaneous  Emission  of  X-Radiation",  Phys.  Rev.  A.,  Vol.  10 
1974,  p.  A903.  

9.  J.  G.  Kepros,  E.  M.  Eyring,  F.  W.  Cagle,  Jr.,  "Experimental 

Evidence  of  an  X-Ray  Laser",  Proc.  Nat.  Acad.  Sci.  USA. Vol.  69.  1972. 

p.  1744.  — — 

10.  R.  A.  Andrews,  "X-Ray  Lasers-Current  Thinking",  NRL  Memorandum 
Report  2677,  October,  1973. 

11.  R.  C.  Elton,  "Analysis  of  X-Ray  Laser  Approaches:  2.  Quasi- 
stationary Inversion  on  Kot-Inncrshel  1 Transitions",  NRL  Memorandum 
Report  2906,  October,  1974. 
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There  are  intriguing  applications  for  an  x-ray  laser.  Reference  10 
lists  and  discusses  a number  of  these  ranging  from  basic  physics  research 
to  medical  studies  but  concludes  that,  at  least  for  early  devices,  x-ray 
lasers  will  have  their  greatest  impact  on  material  science  and  related 
fields.  Some  applications,  however,  would  require  the  beam  to  propagate 
a significant  distance  through  air. 

Unfortunately,  air  is  not  transparent  to  soft  x-ray  photons  as  it 
is  for  optical  photons.  In  fact,  the  range  of  individual  10-keV  x-rays 
in  air  is  only  about  2 meters  and  at  1 keV  the  range  is  smaller  by  a 
factor  of  1000.  However,  these  ranges  are  appropriate  only  for  x-rays 
which  have  a negligible  effect  on  the  medium  through  which  they  propagate. 
If  the  diameter  of  the  beam  is  sufficiently  small  and  the  intensity 
sufficiently  high,  then  the  x-rays  will  have  a drastic  effect  on  the  air. 
In  fact,  there  are  several  effects  produced  which  will  extend  the  range. 
This  study  is  being  made  to  quantify  these  effects  and  determine  x-ray 
laser  parameters  which  maximize  the  range. 


II.  ELEMENTARY  CONSIDERATIONS 

When  a narrow  collimated  beam  of  soft  x-rays  propagates  through 
matter  it  is  attenuated  exponentially  by  absorption  and  scattering 
processes.  The  former  is  due  mainly  to  the  photoelectric  effect  which 
results  in  the  ejection  of  an  electron  from  the  absorbing  atom,  usually 
from  an  inner  shell.  The  atom  then  promptly  fills  this  vacancy  with 
an  electron  from  an  outer  shell  and  the  excess  energy  is  emittec  either 
as  an  x-ray  or,  as  is  usually  the  case  with  atmospheric  gases,  in  the 
kinetic  energy  of  an  Auger  electron.  The  resulting  ion,  which  in  the 
case  of  an  Auger  transition  is  twice  ionized,  may  repeat  the  process  and 
thereby  become  four  times  ionized.  For  nitrogen  (Z=7) , further  photo- 
electric absorption  is  accompanied  by  an  x-ray  rather  than  an  Auger 
emission  so  that  the  successive  absorption  by  an  atom  of  five  x-rays 
can  completely  strip  the  atom  of  electrons.  Further  photoelectric 
absorption  by  this  ion  does  not  occur  until  it  recombines  with  a free 
electron,  and  this  occurs  on  a time  scale  greater  than  probable  pulse 
widths . 

This  stripping  process  is  therefore  one  means  for  reducing  the 
amount  of  photoelectric  absorption,  which  is  the  most  important 
absorptive  interaction  for  soft  x-rays.  Furthermore,  the  amount  of 
energy  expended  by  the  pulse  to  produce  this  stripping  can  be  made  small 
by  making  the  pulse  diameter  small.  This  not  only  reduces  the  number 
of  air  molecules  in  the  path  of  the  pulse,  but  also  reduces  the  distance 
that  an  ion  must  travel  as  a result  of  thermal  motion  in  order  to 
leave  the  path  of  the  pulse.  This  thermal  motion  is  important  because 
the  heating  is  confined  to  a small  region  surrounding  the  beam  and  is 
greatest  within  the  beam.  This  results  in  a large  radial  pressure 
gradient  which  rapidly  removes  particles  from  the  path  of  the  pulse 
because  its  radius  is  very  small.  If  the  x-ray  pulse  has  the  optimum 
shape  in  time,  this  removal  of  absorbing  atoms  can  occur  before  the 
arrival  of  the  main  part  of  the  pulse,  so  that  very  few  atoms  must  be 
fully  stripped  to  eliminate  photoelectric  absorption. 
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When  an  energetic  photoelectron  is  produced  it  starts  away  from 
the  atom  in  a direction  which  would  quickly  remove  it  from  the  path  of 
the  pulse.  The  actual  path  of  the  electron,  however,  is  a tortuous  one 
because  it  can  suffer  large  energy  and  momentum  changes  with  each  collision. 
Even  so,  if  the  air  is  relatively  undisturbed  then  the  range  of  these 
electrons  is  considerably  greater  than  the  pulse  radius.  As  the  air 
then  becomes  slightly  ionized  there  will  be  an  electric  field  produced 
as  a result  of  the  charge  distribution  which  can  roughly  be  described 
as  a small-diameter  cylinder  of  high  density  positive  charge  surrounded 
by  a large-diameter  cylinder  of  low-density  negative  charge.  The  force 
exerted  by  this  field  on  the  energetic  photoelectrons  opposes  their 
outward  motion.  A force  is  also  exerted  on  the  positive  ions  driving 
them  radially  outward.  Tne  resulting  motion  is  a damped  plasma 
oscillation  which  quickly  absorbs  the  excess  kinetic  energy  of  any 
charged  particles,  but  also  raises  the  temperature  and  thereby  leads 
to  a reduction  of  the  number  of  absorbing  atoms  that  the  x-ray  pulse 
encounters,  because  the  net  flow  of  particles  is  radially  outward. 

There  are  thus  at  least  three  processes  which  serve  to  reduce  the  amount 
of  photoelectric  absorption,  which  normally  limits  the  range  of  10-keV 
x-rays,  for  example,  to  about  two  meters. 

But,  even  if  photoelectric  absorption  can  be  neglected  completely, 
the  range  of  x-rays  is  still  limited  to  about  forty  meters  by  the  Compton 
effect,  which  is  scattering  of  x-rays  by  free  or  loosely-bound  electrons. 
However,  the  processes  which  reduce  the  amount  of  photoelectric 
absorption  also  affect  the  amount  of  Compton  absorption.  The  stripping 
of  the  electrons  from  the  atoms  actually  increases  slightly  the  amount 
of  Compton  absorption.  On  the  other  hand,  if  the  net  hydrodynamic  and 
electrodynamic  motion  of  the  positive  ions  is  radially  outward  then  the 
electrons  can't  be  far  behind  because  they  are  much  lighter  and  can  be 
easily  pulled  along  as  the  plasma  tries  to  restore  charge  neutrality. 

There  is  thus  a conceivable  means  for  reducing  both  photoelectric  and 
Compton  losses. 

What,  then,  ultimately  limits  the  range  of  intense  x-ray  pulses? 

As  it  turns  out,  the  pulse  diameter  must  be  quite  small  in  order  to 
reduce  photoelectric  and  Compton  absorption  to  insignificance.  Just  how 
small  it  must  be  is  determined  by  how  efficient  the  hydrodynamic  and 
electrodynamic  motion  is  at  removing  absorbers  from  the  pulse  path. 

The  accuracy  of  such  a prediction  depends  strongly  on  how  well  these 
processes  are  described  mathematically  and  computationally.  Although 
the  formulation  of  a solution  to  these  problems  is  far  from  complete  it 
is  already  apparent  that  the  pulse  diameter  cannot  be  greater  than 
about  10  meters.  But  even  if  the  diameter  is  this  large,  the  inherent 
divergence  of  the  pulse  is  enough  to  reduce  the  intensity  of  a 10-keV 
x-ray  pulse  by  a factor  of  4 as  it  propagates  100  meters.  This  increase 
of  the  diameter  of  the  pulse  reduces  the  rate  at  which  absorbers  move 
out  of  the  path  of  the  pulse,  and  eventually  vhe  diameter  gets  so 
large  and  the  intensity  so  low  that  the  remainder  of  the  range  is  equal 
to  the  normal  range  of  these  x-rays  in  air,  which  is  less  than  two  meters. 
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thi^\°ther  hand’  thi.  ma/  n0t  be  the  ultimate  limit  at  all  if  some- 
rerr«iS?SPf1S>  °T  ^ be  made  t0  happen»  which  focusses  the  x-ray  beam. 

are  St  L IT  "0t  but-  “s  "e  “in  ««>  the” 

SL?  fn"? In -?e  ^ver«ence  angle  which  must  be  overcome  is  only 
'_j:  n , radiana  Because  of  motion  of  ions  and  electrons 

nadxally  outward,  the  late  part  of  the  x-ray  pulse  will  -see'-  mostly 
empty  space,  or  a vacuum,  which  has  a refractive  index  of  unity.  The 

^efractJon^hvf8'^-  late  Part  °f  the  pulse  wil1  have  an  ^dex  of 
rnt  i J!  h h>  1S  eqUal  t0  that  for  normal  air,  can  produce 

°Jce  ZllT^rL  'nCident  on  a sharP  interface  at  the  diver- 

' le;  fhe  complete  description  of  the  conditions  required  for 

noi  ve^  ^en6am  S®}f'focussing  is  a truly  ambitious  undertaking  and  has 
rffecriust?^11^ ta!lvel>\  Sloped,  but  the  possibility  of  such  an 
non-diverging  x^raj  beams!  gatl°n  °f  ^ pr°pagatjon  characteristics  of 


HI. 


NUMERICAL  TECHNIQUE 


4nrl  fWeCaUcc°fithe  IelatiVely  Sh0rt  lifetimes  °f  x-ray  emitting  species 
eaHv  L f i 7 °f  maklng  mirrors  for  x-rays,  it  is  probable  that 
, • ,y-g  neratlon  x-ray  lasers  will  operate  ir.  a pulsed  mode.  Furthermore 
higher  powers  are  possible  in  this  mode  and  we  anticipate  that  these 

rin!npvn^wded  t0  S1^flcantly  extend  the  range.  We  therefore  wish  to 
v t he.waSSaufc  °f  a cylindrical,  narrow-diameter,  finite-length 

",  y pulse  through  a medium.  Since  air  is  mostly  nitrogen  and  oxygen 

The  air  hv°miC  pr°perties  of  these  are  very  similar,  we  will  approximate 

is  divided  S nitrogen  Sas-  For  the  numerical  treatment,  the  pulse 
_ divicied  into  a number  of  sub-pulses,  each  of  which  has  only  a minor 

ettect  on  the  propagation  medium.  The  problem  then  is  to  determine  the 
amount  of  attenuation  suffered  by  each  sub-pulse.  This  attenua^on  lill 
y "J*  °nly  Wlth  the  distance  along  the  length  of  the  pulse  but  also 
with  the  distance  radially  from  the  axis  of  symmetry.  Furthermore  all 
of  these  vary  with  the  distance  the  pulse  has  propagated.  ’ 

nrnmSl"“  the  changes  in  the  attenuation  are  caused  by  changes  in  the 
propagation  medium  it  is  obvious  that  the  medium  must  be  described  by 
some  sort  of  computational  grid.  The  crucial  question  is  whether  this 
1 - VT  °r  three-dimensional,  for  if  the  computer  must 
simultaneously  keep  track  of  the  condition  of  the  medium  along  the 
entire  iength  of  the  pulse  the  memory  and  computation  time  required 
I°U  U 6 VCry  great!  Fortunately  a two-dimensional  grid  ..s  sufficient 
tor  the  x-ray  energies  and  pulse  diameters  with  which  we  are  concerned, 
his  is  true  because  mass  motion  parallel  to  the  direction  of  propagation 

the  JhlS  follows  from  the  fact  that  the  temperature  along 

the.  propagation  direction  changes  more  slowly  with  distance  than  that 

through dlrecti°a'  Furthermore,  the  area  of  the  surface 
irough  which  these  mass  flows  occur  is  much  smaller  for  flows  parallel 

x-ravepSls^intensi?vreCtf01i1,  FdgUre  1 ^ shows  an  assumed  longitudinal 
ra.  pulse  intensity  profile  and  its  approximate  representation  as  a 
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series  of  square  sub-pulses.  Similary,  Figure  1(b)  shows  a transverse 
intensity  profile  which  extends  out  to  some  maximum  radius  rm.  The 
actual  shape  of  the  region  of  space  over  which  the  intensity  is  constant 
is  a long,  thin  cylindrical  shell.  The  entire  pulse  then  looks  like  a 
long,  thin  sliver  and  we  consider  its  propagation  successively  through 
short  intervals  as  shown  in  Figure  1(c). 

We  will  consider  the  propagation  of  the  entire  pulse,  one  small 
section  at  a time,  through  an  interval  of  length  l.  With  the  passage 
of  each  of  these  sub-pulses  a number  of  fundamental  interactions  will 
be  produced  which  alter  the  condition  of  the  medium  contained  in  this 
cylinder.  This  cylinder  is  divided  into  many  cylindrical  shells  of 
equal  thicknesses  Ar  and  outer  radii  rj.  During  a subpulse  the  motion 
of  matter  between  these  cylindrical  shells  is  calculated  from  their 
thermodynamic  properties  assuming  the  material  is  an  ideal  gas.  The 
largest  r^  is  chosen  so  that  the  disturbance  in  the  medium  reaches  the 
outermost  cylindrical  shell  only  after  the  last  sub-pulse  has  passed. 

The  length  of  the  propagation  interval,  l,  is  chosen  so  that  the  overall 
pulse  shape  does  not  change  greatly  from  one  interval  to  the  next. 
Similarly,  the  size  of  the  subpulses  is  chosen  small  enough  so  that  the 
medium  does  not  change  greatly  during  a single  subpulse.  These 
subpulses  must  be  smallest  where  the  absorption  is  greatest  along  the 
pulse,  and  the  relative  position  of  this  maximum-absorption  region  within 
the  pulse  moves  as  the  pulse  propagates.  When  calculations  are  made 
away  from  this  region  several  subpulses  are  combined  to  reduce  computation 
time. 


IV.  MATHEMATICAL  FORMULATION 


A.  Photoelectric  Effect 


The  photoelectric  effect  is  responsible  for  the  severe  attenuation 
of  soft  x-ray  beams  in  matter.  In  this  process  the  x-ray  energy  is  trans- 
ferred to  a bound  electron,  forcibly  ejecting  it  from  the  atom  or  ion. 

The  resulting  ion  may  subsequently  absorb  more  x-rays  by  the  same  process 
until  all  of  its  electrons  have  been  removed.  Photoabsorption  may,  of 
course,  occur  for  electrons  in  the  K-shell  or  in  the  L-shell.  For  (nitrogen) 
neutral  atoms,  photoabsorption  by  K-shell  electrons  is  about  twenty  times 
more  probable.  Since  is  orbitals  are  not  very  sensitive  to  the  occupancy 
of  2S  or  2P  orbitals,  the  cross  section  for  K-shell  photoabsorption  can 
be  assumed  to  be  constant  for  atoms  and  ions  except  when  the  K-shell 
contains  less  than  two  electrons12.  The  L-shell  cross  sections,  however. 


12.  D.  W„  Missavage  and  S.  T.  Manson,  "Photoionization  of  Positive 
Ions  of  Atomic  Oxygen",  Phys.  Letters,  Vol.  38A,  1972,  p.  85. 
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do  vary  with  the  degree  of  ionization  and  are  scaled  from  the  neutral 
atom  values  by  using  Slater's  rules  for  the  screening  potential13.  For 
absorption  by  electrons  in  the  L-shell  the  effective  charge  is 

2eff  - 1 ~ Zk(0.85)  ~ Zl(0.35), 

where  Z is  the  nuclear  charge,  and  Zk  and  Zk  are  the  number  of  electrons 

in  the  K-  and  L-shell,  respectively.  For  the  x-ray  energies  we  consider 

we  will  determine  the  ion  cross  section  per  L-electron  from  the  neutral 

atom  cross  section  assuming  that  the  cross  section  is  proportional  to 

Z 4 
eff  * 

For  each  type  of  absorption  we  will  assume  that  each  subpulse  is 
absorbed  exponentially  as  it  propagates  through  the  distance  l so  that 


where  N0  is  the  initial  number  of  x-ray  quanta  in  the  subpulse,  N is 
the  number  remaining,  y is  the  photoabsorption  cross  section,  and  p 
is  the  density  of  species  in  the  medium  for  which  the  absorption  is 
being  calculated.  We  further  assume  that  the  medium  is  pure  nitrogen  gas 
and  that  the  molecular  bond  is  broken  with  the  first  photoabsorption,, 

Thus  for  each  element  of  the  cylindrical  grid  we  will  keep  track  of  the 
densities  of  nine  different  species:  free  electrons,  neutral  atoms, 

and  ions  with  charges  from  1+  to  7+.  The  cross  sections  for  neutral 
atoms  are  taken  from  analytical  approximations  to  experimental  values14 
which  are  efficient  for  computer  use. 

The  number  of  photoelectrons  produced  in  a pathlength  i is  thus 
(N0-N)  regardless  of  which  component  is  responsible  for  the  absorption. 

The  energies  of  these  electrons  will,  however,  depend  upon  the  degree 
of  ionization  of  the  absorber  since  the  electron  energy  is 

E = hv  - B, 

where  hv  is  the  x-ray  photon  energy  and  B is  the  binding  energy  of  the 
electron  and  varies  with  the  shell  and  number  of  electrons  on  the  ion. 

For  low  energy  such  that  S = v/c  for  the  ejected  electron  is  small,  theory 
and  experiment  are  in  agreement  in  giving  a distribution  in  solid  angle 
of  the  direction  of  ejection  of  the  photoelectrons  proportional  to 


13.  M.  D.  Torrey,  N.  A.  Zessoules,  and  W.  D.  Lanning,  "Numerical 
Simulation  Techniques  for  Nonequilibrium  Effects  of  High-Aititude 
Nuclear  Detonations",  Technical  Operations  Research,  Report  No.  TO-B 
66-29,  May,  1966. 

14.  F.  Biggs  and  R.  Lighthill,  "Analytical  Approximations  for  X-Ray 
Cross  Sections  II",  Sandia  Laboratories,  Report  No.  SC-RR-71  0507, 
December,  1971. 


± 


cos2')',  Where  4-  is  the  angle  between  the  direction  of  ejection  and  the 
e ectric  vector  of  the  photon  , which  is  perpendicular  to  the  propaga- 
tion direction.  We  will  simply  assume  that  all  photoelectrons  are 
emitted  in  the  most  probable  direction,  which  is  radially  outward  from 
the  x-ray  path.  We  also  assume  that  each  ion  is  in  its  ground  state 
before  the  absorption  and  returns  to  the  ground  state  of  the  new  ion 
promptly  after  the  absorption.  Since  the  fluorescence  yield  of  nitrogen 
is  only  about  .00619,  this  is  tantamount  to  assuming  that  an  Auger 
emission  promptly  follows  the  K-shell  absorption  by  ions  which  have  four 
or  more  electrons  and  that  a low-energy  x-ray  is  promptly  emitted  after 
a three-electron  ion  absorbs  an  x-ray  in  the  K-shell.  The  energy  of  Auger 
electrons  is,  to  a good  approximation, 

E » W(K)  - 2 W(L) , 

iS,  the  ener8y  of  the  ion  in  which  a K-electron  is  missing 
and  W(L)  is  the  energy  if  an  L-electron  is  missing17.  The  energy  of  the 
uger  electrons  is  much  less  than  that  of  the  photoelectrons  and  results 
in  heating  much  closer  to  the  propagation  path. 

h*  Energetic  Electron  Energy  Loss 


When  an  atom  or  ion  absorbs  an  x-ray  via  the  photoelectric  effect 
the  x-ray  disappears  completely  and  its  energy  is  transferred  to  a 
?ound  electron.  Part  of  this  energy  is  transferred  to  the  ion  as  its 
iond  with  the  electron  is  broken.  The  remaining  energy  is  carried  by  the 
photoelectron  as  kinetic  energy  and  is  deposited  in  the  medium  some 
distance  trom  where  the  primary  interaction  occurred.  The  size  of  the 
region  where  the  photoelectron  loses  its  energy  and  the  manner  in  which 
it  does  so  depends  upon  the  condition  of  the  medium.  When  the  medium  is 
in  its  initial  state  this  energy  loss  proceeds  by  a well -understood 
mechanism  and  has  been  extensively  studied  experimentally.  However,  as 
the  air  becomes  ionized  and  the  charge  separation  begins  to  produce  an 
electric  field,  the  energy-loss  mechanism  is  not  so  clear  and  experimental 
data  is  not  available.  On  the  other  hand,  some  simplifying  assumptions 
can  be  made  which  lead  to  a representation  of  the  problem  which  is  amenable 
to  a computer  calculation. 


..l5,  A’_  H*  ComPton  and  S.  K.  Allison,  X-Rays  in  Theory  and  Experiment 

Van  Nostrand,  Princeton,  N.J.,  1935,  pp.  564-582.  

r n ?am5ynik*  B’  baseman,  R.  W.  Fink,  H.  U.  Freund,  H.  Mark, 

‘ * wift,  R.  E.  Price,  and  P.  V.  Rao,  "X-Ray  Fluorescence  Yields, 

Auger,  and  Coster-Krorig  Transition  Probabilities",  Rev.  Mod  Phvs 
Vol.  44,  1972,  p.  716.  ~ 

1 7.  E.  U.  Condon,  "X-Rays",  in  Handbook  of  Physics  edited  by  E.U. 
Condon  and  H.  Odishaw,  McGraw-Hill,  New  York,  1958,  p.  7-124. 
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In  a relatively  undisturbed  medium  energetic  electrons  lose  energy 
primarily  via  ionizing  collisions  with  neutral  atoms.  The  energy  loss 
per  unit  path  length  for  this  process  is  given  by18 


where  r0  = ^ ® classical  electron  radius 

4nro2  B 1.00  x 10  2i*  cm2/electron  * 1.00  barn/electron 
MeC2  = 0.51  MeV  « electron  rest  energy 
62  = (v/c)2  = 1-  [ (T/M0c2)  + 1]~2 
NZ  = bound  electron  density 
I = geometric  mean  ionization  and  excitation  potential 
The  energy  I can  be  written 


Z In  I 


In  A . , 
n,  1 


(2) 


where  f„  j is  the  sum  of  the  oscillator  strengths  for  all  optical  transi- 
tions of  the  electron  in  the  n,l  shell  and  is  close  to  unity  while  An  ^ 
is  the  mean  excitation  energy  of  the  n,l  shell  and  can  be  put  equal  ’ 
to  the  ionization  potential  with  sufficient  accuracy.  Theoretical  values 
of  I,  calculated  with  hydrogen-like  wave  functions,  are  of  little  practical 
value.  But  empirical  values  of  I determined  for  each  value  of  Z can  be 
used  to  accurately  predict  the  energy  loss  in  unionized  matter19. 

However,  in  the  problem  at  hand,  we  must  be  concerned  also  with 
the  further  ionization  of  ions.  The  value  of  I will  certainly  be  different 
for  these,  no  experimental  measurements  are  available,  and  no  rigorous 
prescription  exists  for  determining  I.  But  rather  than  neglect  the 
variation  of  I with  the  state  of  ionization,  we  shall  employ  equation 
(2)  taking  fn>i  = 1,  An>i  = ionization  potential,  and  summing  only  over 
those  electronic  states  occupied  when  the  ion  is  in  its  ground  state. 


18.  R.  D.  Evans,  The  Atomic  Nucleus.  McGraw-Hill.  New  York.  1955. 
p.  532.  — " 


I9.  Ibid,  p.  580. 
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for  nitrogen,  this  results  in  values  for  I ranging  from  670  eV  for  the 
b+  ion  to  92  eV  for  neutral  atoms.  This  can  be  compared  to  the  experi- 
mental value  for  neutral  nitrogen  which  is  80.1  eV2°. 

As  the  medium  becomes  ionized  the  free  electron  density  increases  and 
energetic  electrons  will  lose  energy  in  collisions  with  these.  Follow- 
ing the  approach  taken  by  Richards  , we  will  write  the  energy  loss  for 
fast  electrons  in  a plasma  as 


dE 


dT 


6'4  * Mas 


\(0 


ion 


(3) 


where  the  first  term  accounts  for  the  energy  lost  by  electrons  to  other 
electrons;  ne  is  the  free  electron  density  and  $c  is  the  Coulomb  cross 
section,  given  approximately  by  -rr(e2/E) 2 . The  second  term  is  the  sum 
of  the  energy  losses  due  to  ionizing  collisions  with  each  ionic  species 
computed  from  equation  (1).  The  last  term  accounts  for  short-range 
collisions  and  for  excitations  of  the  ions  and  is  assumed  to  be  equal 
to  the  middle  term,  as  it  is  for  neutral  gases. 


Another  phenomenon  with  which  we  must  be  concerned  is  the  conversion 
of  the  kinetic  energy  of  photoelectrons  to  potential  energy  in  an 
electric  field.  Although  this  is  an  important  effect,  its  rigorous 
inclusion  would  complicate  and  greatly  increase  the  computation  time  of 
the  problem.  The  effect  of  this  phenomenon  will  be  to  confine  the 
region  of  photoelectron  energy  loss  to  the  immediate  vicinity  of  the 
plasma-which  in  turn  will  be  near  the  beam.  In  order  to  account  for 
this  at  least  in  an  approximate  way,  we  will  assume  that  photoelectrons 
lose  their  energy  in  the  cylindrical  shell  in  which  they  are  produced. 

C.  Compton  Effect 


The  Compton  effect  produces  both  scattering  and  absorption.  That 
is,  part  oL  the  energy  of  the  affected  photons  is  transferred  to  the 
responsible  electron  and  adds  to  the  thermal  energy.  However,  most  of 
the  energy  is  scattered  and  the  photon  has  no  further  appreciable  effect 
on  the  medium.  This  loss  mechanism  therefore  can  use  up  large  amounts 
of  the  x-ray  pulse  energy  and  limit  the  range  to  about  40  meters. 


20*  Kai  Siegbham  Alpha-,  Beta-,  and  Gamma-Ray  Spectroscopy.  North- 
Holland,  Amsterdam,  1966,  p.  26.  

21.  P.  I.  Richards,  "Summary  Report  on  Investigation  of  Radiation 
and  Chemical  Calculations",  Technical  Operations  Research,  Report  No.  TO-B 
62-24,  May,  1962.  r 
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Compto^ran ee  CCtrlC  k^86  i$  C0nsic3erably  less  than  the 

’•  H h medlum  Can  be  Sreatly  disturbed  before  Compton 

the  elertr0meHimPOrta^t’  the  nature  of  tbe  disturbance  is  to  reduce 
can  be  redTced  Y ° Vicinlt>'  of  the  then  Compton  losses 

The  average  Compton  collision  cross  section  is  given  by 

8tt  o 

e°  ~ T ro"  Cl'^a  + 5,2a2  - 13.3a3+...)  (4) 

for  small  a - _£vo,  where  hv0  is  the  x-ray  energy  and  M0c2  is  the 

It ■e!t2n/BSi  enerSy22-  ^is  cross  section  is  the  probability  for  removal 

cont,?n  mnoner°emiea.COllr5ed  »hile  Pass-8  ‘hSough  an  absorber 
is  giVen  bv°3  electron/cm  • rh£;  average  Compton  absorption  cross  section 


% 

0 = 
e a 


(a-4.2a2  + 14.7a3-...) 


The  average  energy  removed  from  an  incident  beam  by  Compton  collisions  is 
measured  bv  eo  while  the  average  energy  absorbed  by  the  electrons  is 
measured  by  The  average  energy  of  a Compton  electron  is  then 


This  amount  of  energy  will  be  added  to  the  total  thermal  energy  for  each 
Compton  electron  and  will  thus  appear  directly  as  heat.  Other  explicit 
processes  which  contribute  directly  to  the  thermal  energy  include  recoil 
of  an  ion  with  photoelectron  emission,  energy  transferred  in  photoelectron- 
thermal  electron  collisions,  and  the  energy  of  photoelectrons  which  is 
less  than  the  lowest  ionization  potential  but  greater  than  the  thermal 
energy. 

0.  Mass  Motion 

pie  absorption  by  a medium  of  a narrow  x-ray  beam  quickly  produces 
a cylinder  which  has  a large  radial  temperature  gradient  and  a negligible 
axial  temperature  gradient.  The  resulting  mass  motion  is  of  utmost 
importance  since  it  removes  absorbing  particles  from  the  path  of  the  beam, 
ror  small  beam  diameters  this  can  occur  rapidly  enough  for  the  resulting 
range  to  be  considerably  greater  than  it  would  be  otherwise. 


22.  R,  D.  Evans,  op.  cit.,  p.  684. 

23.  Ibid,  p.  688. 
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In  the  limit  as  one  goes  to  smaller  beam  diameters  and  shorter 
pulse  widths  this  mass  motion  changes  from  a hydrodynamic  flow  to  an 
effusion  flow  because  the  heated  ions  can  move  distances  significant 
compared  to  the  beam  diameter  before  they  equilibrate  with  the  surround- 
ing medium.  One  may  then  assume  that  the  heated  region  is  a radially 
expanding  cylinder  with  its  escaping  particles  effusing  into  the 
surrounding  medium  and  gradually  raising  its  temperature.  Calculations 
have  been  made  with  this  approximation  and  the  results  are  given  in 
the  next  section. 

A more  exact  representation  of  the  mass  motion  can  be  made  by 
including  the  physics  of  fluid  flows.  Following  a method  used  to 
describe  explosive  hydrodynamics24,  we  employ  a finite-difference 
approximation  to  the  differential  equations  which  expresses  the  basic 
physical  limitations  on  the  flow  without  introducing  spurious,  unstable 
processes  through  the  numerical  approximations. 

The  fundamental  physical  laws  of  fluid  flow  are  conservation  of 
mass,  momentum,  and  energy.  The  state  of  the  fluid  can  be  specified 
with  the  densities  of  each  of  these  quantities.  The  conservation  laws 
are  imposed  by  requiring  that  flows  out  of  one  cell  of  the  computational 
gud  must  reappear  in  another.  This  is  most  easily  ensured  by  focussing 
attention  on  the  flows  across  the  faces  between  grid  cells.  Computing 
these  interfacial  flows  for  the  rpass  essentially  evaluates  V • (pv) , 
where  p is  the  mass  density  and  v is  the  fluid  flow  velocity.  Crediting 

and  debiting  these  flows,  for  a brief  step  forward  in  time,  essentially 
sets  ' 


ff-=-V-(pv),  (7) 

which  expresses  the  law  of  conservation  of  mass. 

A coi responding  treatment  for  the  energy  expresses 

ff“  -V.  [d>p)v] 

where  p is  the  pressure  and  E = e + pv2, 

~ 

e is  the  thermal  energy  density,  and  E^is  the  total  energy  density. 
Gravity  is  neglected  and  the  term  V« (pv)  accounts  for  the  work  done 
by  compression  and  by  streaming  against  the  local  pressure  gradient. 
Ihe  net  result  is  that  one  must  evaluate  interfacial  flows  of  the 
quantity  E+p. 


(8) 

(9) 


2^.  P.  I.  Richards,  "The  Tech/Ops  Two-Dimensional  Time-Dependent 
Hydrodynamics  Code",  Technical  Operations  Research,  Report  No  TO-B 
62-22,  April,  1962.  v 
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is 


Since  the  axial 
expressed  as 


motion  is  ignored,  the  momentum  conservation 


law 


at 


* ~^*  [(pv  )v]  - |E 
r Br 


A straightforward  time  integration  of  the  form 


XCt+At)  = y (t)  + ijl  At 


(10) 


(11) 


t]?  Chf1ReS  in  the  mass*  enersy*  and  momentum  densities 
he  of  the  time  steps  At  corresponds  to  the  widths  of  x-ray  sub- 

At  isSsmanmrnmn^e  7+U  7°^  *°  ^ the  maxl,num  sWl  velocity  times 
At  is  small  compared  to  the  computational  grid  dimension, 

E.  Other  Processes 

Multiple  ionisation  may  occur  following  K-shell  photoelectric 
absorption  due  to  a process  called  electron  shake-off.  This  is  the 
phenomenon  by  which  an  electron  in  a given  orbital  is  excited  into  a 
new  orbital  or  into  the  continuum  as  the  result  of  a sudden  change  in 

m ofeJhJal  P?tentlal°  Since  the  Potential  seen  by  the  electron  is  made 
‘wt  nuclear  charge  minus  the  shielding  of  the  other  electrons, 

nuifiir  ,shdke“off  ma^  be  initiated  either  by  a sudden  change  in  the 
n clear  charge,  as  m 0-decay,  or  a sudden  change  in  the  electron 
configuration,  as  in  photoionization.  The  most  successful  treatment  of 

aid  , ;°f5  ,haS  been  throu8h  the  use  of  the  sudden  approximation 

ia  i f shake-off  probabilities  following  0-decay  have  been 

f f.  f°r  many  clements  • However , calculations  have  not  been 

f i -I  ?r  b particular  case  in  which  we  are  interested:  shake-off 

his  heiiS-t,  H^°a°nizatl°n  ln2pip708eT1,  The  nearby  element  neon,  however, 
has  been  studied  extensively2* » 27  and  experimental  data  is  also 


-5.  T„  A.  Carlson,  C.  W.  Nestor,  Jr.,  T. 
B«  Malik,  "Calculation  of  Electron  Shake-off 
Z-2  to  92  with  the  use  of  Self-Consistent-Field 
IVA-  Rev.,  Vol.  169,  1968,  p.  27c 


C.  Tucker,  and 
for  Elements  from 
Wave  Functions", 


n m26'c,  T’  A'  Carlson  and  M*  0.  Krause,  "Experimental  Evidence 
^oubie^Electron  Emission  in  an  Auger  Process”,  Phys.  Rev.  Letters 
vol.  14,  1965 , p„  390.  —t — 


for 


27,  M.  0.  Krause,  M.  L.  Vestal,  W.  H. 
Readjustment  of  the  Neon  Atom  Ionized  in 
Phys,  Rev.  A,  Vol.  133,  1964,  p.  A385. 


Johnston,  and  T.  A.  Carlson 
the  K Shell  by  X-Rays", 
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avaiiabie  for  the  dependence  of  the  shake-off  probability  on  the 
photoelectron  energy  . The  combined  probability  for  the  shake  off 

e«L°rir,rbo^ir«%ln  additi0"  “ thaphoS^'Sro^aSt^f 

electron  is  about  25. 4t  for  a photoelectron  energy  of  7.5  keV. 

„ ibbikebion  potential  of  an  ion  in  a plasma  is  not  the  same 

“,,$9  “It^hroosU^V111  redU“d  beC“USe  of  PlasI"a  P°>arira- 

times  the  ele«ro“J  charge  “the0"  "Uh  ionic  ct!a^=  <=<iuai  to  R 
chareed  T8'  the  average  potential  from  all  other 

solufion  of  P^ss  " 6 sy?tem  iS  ob,tai"ed  from  an  approximate 

.solution  of  Poisson’s  equation  and  is30 

J’CR)  = 9 [2  (Z*+l)]-l  ( [3(Z*+1)  KCR)  + I]2/3  - 1},  fl2) 

where  6 is  the  temperature  in  energy  units.  The  other  Quantities 

ZsitTS  " f“8ti°n  02)  are  de«"=d  d"  «™s  o?  tSe  av“aje  electron 
density  Ne  ami  the  average  density  NR  of  ion  species  with  charge  R, 


Z ■ ~ - £r  NR  «2/Ne  (13) 

K(R)  = Re2/ (D0)  C14) 

D"2  = 4tt  Ng  e2 (1  + Z*) /Q  (15) 

ZZt  hTmade”^  “^KgTcR?  £ ” ° * ^^amically 

d(R)  • (R/zj  J-  fT).  (16) 

lower  than  tha/of^^oUted  ioH  by  tlirfmouni"  th<!  PlaS”a  ^ the" 


f™  ?;h.n  r'  CarlS°"  8nd  ”•  0<  Krausa-  "Electron  Shake-off  Resulting 

W5SWJS; " ; •'  “™ 

NLL  and  S0F™ARt,  Report  No.  DASA  2597-1,  April,  1971. 
in  Plasmas",  Ap.  J..  Vol^lJl*  1966^^1  ZO^1™8  °f  l0IUZatl0n  Potentials 
Constants  of 
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where 


Mr  * (R+l)  J (Z) , 
T=  n"1  n 

e 


U = N"1  l Nu  R2 
R R 
and 


Electron-ion  recombination  and  the  resulting  radiation  can  alter 
the  absorption  characteristics  for  long  width  pulses.  This  process 
leads  to  the  reappearance  of  photoelectric  absorbers  on  a time  scale 
determined  by  the  recombination  coefficient,  theoretical  expressions 
for  which  are  given  by  Spitzer32.  These  theoretical  values,  however, 
are  usually  many  orders  of  magnitude  smaller  than  measured  values  in 
laboratory  plasmas.  This  apparent  discrepancy  may  be  attributed  to  a 
different  type  of  recombination  process33*34  in  which  a molecular  ion 
captures  an  electron  and  dissociates  into  two  neutral  atoms.  Quantum- 
mechanical  calculations  predict  an  electron-ion  recombination  coefficient 
of  the  order  of  10"12  cnr/sec35.  This  indicates  that  recombination 
needs  to  be  considered  for  ions  which  remain  exposed  to  the  pulse  for 
periods  greater  than  or  about  100  ns.  However,  for  the  intensities 
and  densities  currently  being  studied,  recombination  can  be  ignored. 

Several  other  phenomena  are  also  not  included  in  the  calculation 
at  this  stage.  Among  these  are  the  effects  of  induced  electric  and 
magnetic  fields,  inverse  Bremsstrahlung  absorption,  beam  divergence,  and 


32.  Lyman  Spitzer,  Jr.,  Physics  of  Fully  Ionized  Gases,  Inter- 
science, 1962,  p„  150. 

33.  D.  R.  Bates,  "Electro*  Recombination  in  Helium,"  Phvs.  Rev., 

Vol.  77,  1950,  p.  718.  * 

34.  D.  R.  Bates,  "Dissociative  Recombination",  Phys.  Rev.,  Vol.  78, 
1950,  p.  492. 

35.  E.  C.  G.  Stueckelberg  and  P.  M.  Morse,  "Computations  of  the 
Effective  Cross  Section  for  the  Recombination  of  Electrons  With 
Hydrogen  lens",  Phys.  Rev.,  Vol.  36,  1930,  p.  16. 
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haie"beenSm^P  whf1Jh°U8h  th6Se  aTG  n0t  included  explicitly,  assumptions 
een  made  which  are  expected  to  produce  results  qualitatively 

simi  ar  to  those  which  would  be  obtained  from  an  exact  treatment.  In 
part iciaar  the  effect  of  an  electric  field  would  be  to  greatly  shorten 
t range  o photoelectrons  ejected  radially  outward  from  the  beam  path. 
Hence  the  assumption  is  made  that  photoelectrons  are  absorbed  in  the  com- 
putational grid  in  which  they  are  produced. 

The  inherent  beam  divergence,  X,  can  be  large  enough  to  considerably 

spread  the  beam  when  the  diameter  is  very  small.  If  this  spread  does 
indeed  occur,  then  the  range  is  greatly  reduced  from  that  expected  for 
parallel  beams.  However,  the  disturbed  medium  has  characteristics  which 
could  lead  to  confinement  of  the  beam  or  even  whole-beam  self- focusing. 

j density  of  the  medium  seen  by  the  latter  part  of  the  pulse  will  be 
reduced  by  a factor  of  MO4.  If  this  can  be  approximated  by  a vacuum 
with  an  index  of  refraction  of  unity,  and  the  surrounding  medium  could  be 
approximated  by  normal  air  with  an  index  of  refraction  of  1-2.33  x 10“9 
tor  10  keV  x-rays,  then  the  grazing  angle  for  total  reflection  at  a 
sharp  boundary  would  be  MO"4.  The  diffraction-limited  divergence  A 

for  this  case  would  be  MO" 5 , indicating  that  x-rays  in  the  latter  part 
of  the  pulse  diverging  at  this  angle  would  be  totally  reflected.  Until 
both  beam  divergence  and  self- focusing  can  be  included  in  a consistent 
manner  we  will,  therefore,  consider  only  parallel  or  non-diverging  beams. 


V.  PRELIMINARY  RESULTS 

Calculations  have  been  made  with  a preliminary  version  of  the 
computer  code  for  a specific  case.  The  laser  parameters  and  predicted 
range  for  this  problem  are  given  in  Table  I.  The  pulse  energy  assumed 
corresponds  to  about  10%  of  the  maximum  energy  which  can  be  stored  in 
a normal  density  solid  cylinder  of  the  specified  diameter  and  one  meter 
m length,  assuming  each  atom  stores  10  keV.  The  range  here  is  defined 
as^the  distance  at  which  the  pulse  intensity  is  attenuated  by  the  factor 
e - or  the  same  attenuation  factor  that  an  x-ray  beam  normally  experi- 
ences m one  mean  free  path.  The  overall  attenuation,  however,  is  not 
an  exponential  function  of  distance  but  is  almost  linear. 

This  result  was  obtained  using  the  effusion  approximation  for  the 
transverse  motion  of  ions  and  electrons.  The  densities  of  four  of 
the  species  as  well  as  the  x-ray  flux  are  shown  in  Figure  2 as  a function 
ot  time  for  the  initial  propagation  interval.  The  x-ray  intensity 
begins  abruptly  five  half-widths  (100  ns)  before  the  arrival  of  the  peak 
and  has  gaussian  shape  and  is  not  visible  on  the  linear  scale  until 
approximately  50  ns  later.  The  computational  grid  in  this  case  consists 
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°!  only  three  zones  and  the  results  nVO 

shown  in  the  figure,  the  density  of  L f 1 10  central  zone.  As 

the  population  of  the  variotK^lorf hinhw!  falls  to  zero,  and 

times.  The  decay  of  the  7+  ionic^ 1 }l^hly-ionized  species  peaks  at  later 

the  path  of  the  beam.  By  the  time'  the1^  ^ ^ t0  maSS  motjon  out  of 
t - 100  ns,  virtually  all  of  til  V peak  X~T'dy  flux  "drives,  at 
and  electrons  have  e'scaped  from  theT  SJriP,ed  and  most  of  the  ions 
amount  of  Compton  absorption  T1  1 hcu*?ath  the"eby  reducing  the 
considerably  greater  than  the  Cometon^  ting  range  CO. 35  km)  is  therefore 
electric  range  (0.002  km)  1 range  (0-04  km^  a"d  the  photo- 

through  t^x-ray  DuSeT*1”  °btJined  in  the  ce"<-al  M region 

warm  slieath  surroundxng  TkeVlT^  ™ ^ W regioa  which  ^ a 
from  region  to  region  is  obtained  motlon  of  the  particles  and  energy 

outrard  at  a rate  ^proportional  "to 

minus  a correction  for  particles  i 7 and  average  velocity 

assumption  is  made  that  the  hivh  nms!  ‘ at  th°  boundary-  'Hie  implicit 
keeps  particles  from  ^iJ^Snd so tZ  gme™tod  in  tbe  M region 
for  long  pulse  widths  Another  1 ■ ? results  are  not  applicable 

ions  move  a considerable  distance  assumPtlon  is  that  the  effusing 

density  wall  docs  nofa^r  aT  h e Zn/^Z"  *° 
results  are  not  applicable  for  larger  pufie  lad  w^Ttl  “ tHc 

as  shown  in  figure  3 ?Ms  Ml?  f,UnCtlon  °f  Propagation  distance, 
the  medium  becomes  transparent1  after  °?C?  Propagation  interval  in 
amount  of  absorption  a certain-  relatively  constant 


Pulse  energy 

Pulse  axial  shape 

Pull  width  at  half  maximum 

Pulse  radial  shape 

Diameter 

Range 


1000  Joules 
Gaussian 
40  ns 

uniform  intensity 
10  u 
0.35  km 


Table  1.  Assumed  Laser  Parameters  and  Predicted  Range. 


fluidtf!w1™nMunnrrcomSt”io„rr«-,Ig  th°-la?s  of  Mrodynamic 
cylindrical  shells  with  thicknesses  of  1 u ld  consibtinS  of  concentric 
this  treatment  is  considerably  greater  than  that  tomPutatlon  tune  using 
approximation.  Results  have  hoen  oi  • an  tbat  usmg  the  effusion 
distribution  in  the  Ss^rop^atio  T* 

distribution  at  two  times  for  the  nterval*  figure  5 shows  this 

Times  t,  and  t2  are  80  ns  and  60  ns'  respe(toel(S  (T  in,  TaMe  ’■ 
of  the  peak  of  the  x-ray  pulse  Th?  J ????  X.  before  the  arrival 

only  a very  small  fraction  the  Dulse  th  ^ corresP°nd  to  passage  of 
interval  and  required  ” T*1  *5°  first  Propagation 

might  be  expected,  computations  ha(e  not  ye)  bL^TT  l°",  As 

distances  so  accurate  predictions  of  the  range  from  this  mert^faS^of*10" 
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TEMPERATURE  (*V)  PULSE  ENERGY  {Joules} 


wi  " "UP 
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Figure  3.  Remaining  pulse  energy  versus  propagation  distance  in  the 
effusion  approximation. 


0 100  200 

TIME  (ns) 


Figure  4.  Temperature  versus  time  in  the  effusion  approximation  within 
the  beam  (region  M)  and  surrounding  the  beam  (region  W) . 
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thre=aord;T,HrVer’-,S^Ce,thC  ,°Veran  densiU  is  redu««  by  more  than 

exnect  that  thf  ”f  °nf  before  tho  F'Ulsc  peak  arrives  one  may 

expe ct  that  the  range  should  be  at  least  1000  times  greater  than  the 
normal  range,  or  greater  than  2 km. 

Figure  6 shows  the  result  of  a similar  calculation  except  that 

1 ximum'nf  ?n  !aken  t0  bc  Gau^im,  with  a full  width  at  half 

maximum  of  10  y The  density  reduction  at  the  center  is  not  as  great 
m tin  s case  and  t lie  range  is  probably  less. 


VI.  SUMMARY 

Several  non-linear  processes  have  been  discussed  which  are  pro- 
duced as  a result  of  the  propagation  of  x-rays  through  air  and  yet 
reduce  the  amount  of  absorption  suffered  by  later  x-rays  which  travel 
the  same  channel.  The  most  important  of  these  processes  is  the  transverse 
hydrodynamic  motion  of  the  gas  particles.  For  typical  x-ray  pulse 
parameters  the  gas  particle  density  can  be  reduced  by  an  order  of 
magnitude  in  a few  nanoseconds.  The  efficiency  of  this  density  reduction 

and  range  extension  depends  upon  the  parameters  which  describe  the 
x-ray  pulse. 

Future  efforts  on  this  study  will  be  aimed  at  including  explicitly 
the  formation  of  electric  fields  and  their  effect  on  the  motion  of  the 
photoelectrons  and  therefore  on  the  size  of  the  region  which  is  heated  by 
,he  beam.  In  addition,  a method  for  describing  diffraction  and  refraction 
. x-rays,  similar  to  methods  applied  to  optical  lasers36*37,  will  be 
included  to  determine  whether  self- focusing  can  be  produced.  Finally, 
calculations  will  be  made  for  a variety  of  x-ray  pulse  parameters  to' 
determine  the  ones  which  optimize  the  range. 


36.  B.R.  Suydam,  "A  Laser  Pulse  Propagation  Code",  Los  Alamos 
Scientific  Laboratory  Report  LA-5607-MS,  1974. 

O)"7.  C.J.  F.liiot  and  D.  13 . Henderson,  "A  New  Technique  for  Nonlinear 
Optical  Problems",  Los  Alamos  Scientific  Laboratory  Report  LA-5639-MS, 
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